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We present the first rigorous derivation of a number of universal relations for a class of models 
with continuously varying indices (among which are interacting planar Ising models, quantum spin 
chains and ID Fermi systems), for which an exact solution is not known, except in a few special 
cases. Most of these formulas were conjectured by Luther and Peschel, KadanofT, Haldane, but only 
checked in the special solvable models; one of them, related to the anisotropic Ashkin- Teller model, 
is novel. 
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It has long been conjectured, mainly by Kadanoff [10- 
12], Luther and Peschel [17] and Haldane [13], that a 
number of universal relations among critical exponents 
and other observables hold in a wide class of models, 
including planar Ising-likc models with quartic inter- 
actions, vertex or Ashkin- Teller models, quantum spin 
chains and ID fermionic systems. Such relations ex- 
press how the universality principle works in models with 
continuously varying indices: the critical exponents are 
model dependent (non-universal) but satisfy model inde- 
pendent formulas, so allowing, for instance, to express all 
the exponents in terms of a single one. The universal re- 
lations have been verified only in certain special exactly 
solvable models, but the conjecture is that they are gen- 
erally valid in a larger class of models, for which an exact 
solution is not available. 

The interest in this kind of universal relations has been 
renewed by recent experiments on materials described 
by models in this class, like quantum spin chain models 
(KCuF3) [16], carbon nanotubes [1], layered structures 
[2] or even ID Bose systems [15]. In such systems the 
critical exponents depend on the extraordinarily com- 
plex and largely unknown microscopic details of the com- 
pounds, but the universal relations allow concrete and 
testable predictions for them in terms of a few measur- 
able parameters. 

Several attempts in the last thirty years have been 
devoted to the proof of the universal relations [19, 20, 
23, 24], by taking as a starting point the formal contin- 
uum limit (identical for all the models considered here), 
where extra Lorentz and Gauge symmetries are verified 
and make it solvable. Of course, lattice effects destroy 
such symmetries and change the exponents; however, this 
problem has never been analyzed. On the other hand, 
not all the relations which are valid in the special solv- 
able models are generically true; this happens, for exam- 
ple, for the exponents involved in the dynamic correla- 
tions [15] and another example will be shown below. It is 
therefore important to determine rigorously, and there- 
fore unambiguously, under which conditions and which 
one among the relations valid in the solvable models are 
generically true. 

Aim of this letter is to report the first rigorous deriva- 
tion of several of such universal relations in a wide class 



of models, including non solvable models; in addition we 
will also prove a relation which is totally new. 

The simplest class of models in the class of universality 
we are considering is coupled Ising models. A configura- 
tion (cr, cr') is the product of two configurations of spins 
cr = {cTx = ±1}xgA and cr' = {cr^ = ±l}xeA- For a 
finite lattice A, the energy i?(cr, cr') is a function of the 
parameters J > 0, J' > and A 

H=-JY. '^-^-+-. - E ^>x+e, - AF(a, CJ') (1) 

i=o,i i=o,i 

where ep and ei are the horizontal and vertical unit bond. 
V{a, cr') is a quartic interaction, short ranged and sym- 
metric in the exchange a ^ a' ; for instance 

^(f^>cr')^ ^ ^ i;(x -y)crxcrx+e,cr^cr^+e, 

j=0,lx,yeA 

with w(x) a short range potential. It is well known that 
several models in Statistical Mechanics can be rewritten 
as coupled Ising models. In particular the Ashkin-Teller 
model [3], a natural generalization of the Ising model 
to four states spins, can be rewritten in the form (1) 
with w(x) = i5x,o- Another example is provided by the 
Eight Vertex model, in which J = J' and V{a, cr') = 

Z]j=0,l ExGA'^x+j(eo+ei)'^x+eoO'x+j(eo+ei)'^x+ei- An CX- 

act solution [3] exists only in the case of the 8V model and 
not for the generic Hamiltonian (1); even in the case of 
the 8V model, the correlations have not been computed 
and only a few indices can be obtained. 

Recently new methods have been introduced in [22] 
and [18] to study 2D statistical mechanics models, which 
can be considered as a perturbation of the Ising model. 
These methods take advantage of the fact that such sys- 
tems can be mapped in systems of interacting fermions 
in d = 1 + 1 dimensions. This mapping was known since 
a long time [21], but only in recent years a great progress 
has been achieved in the evaluation of the Grassmann 
integrals involved in the analysis of the interacting mod- 
els, in the context of Quantum Field Theory and Solid 
State Physics, so that one can take advantage of this new 
technology to get information about statistical mechan- 
ics models. At the moment, when an exact solution is 
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lacking, this is the only way to get rigorous quantita- 
tive information on the properties of such systems. The 
algorithm is based on midtiscale analysis and allows us 
to prove convergence of several thermodynamic functions 
and correlations up to the critical temperature; essential 
ingredients of the analysis are compensations due to the 
anticommutativity of Grassmann variables and asymp- 
totic Ward Identities (WI). 

By using such methods, it has been proved in [18], in 
the case J = J' and A small, that the model is critical 
in the thermodynamic limit at the inverse temperature 
13c = = arctan(\/2- 1)/J-|-0(A); for T near Tc, the 
specific heat behaves as 



]T-Te|-"-l] 



(2) 



with a a continuous non trivial function of A. Moreover, 
if G'^(x — y), £ = ±, are the correlation functions of 
the two quadratic observables = X^j-=o,i i^xCTx+e^ + 
^ Sj=o 1 "'x'^x+ej (which are called energy, if £ = +, and 
crossover, if £ = — , in the AT model, while the names 
are exchanged in the 8V model), in [18] it has been also 
proved that the large distance decay of ^^(x— y) is faster 
than any power of ^^^|x — y|, with correlation length 

^~|r-re|-% asT^Tc 

while at T = Tc, the decay of G^(x — y) is power law: 

G^(x-y)~ |x-y|-2-^ . 

In the Ashkin- Teller model with J ^ J' , it has been 
proved in [9] that there are two critical temperatures, Tc,i 
and Tc 2, such that 



-A-« log[A-2[T - Te,i| • \T - TcaW (3) 



where 2A2 = (T - Tc,if + (T - Tc,2f (the index a in 

(3) is the same as in (2) ). While in the isotropic AT 
the logarithmic singularity of is turned by the inter- 
action in a power law, in the anisotropic AT Cy has still 
a logarithmic singularity; however, Ti^c — l2,c scales with 
a transition index a^T = 1 + 0(A) in the isotropic limit: 



\Tl,c — T2,c\ 



\J-J' 



i\XT 



(4) 



The existence of Xt was overlooked in the literature. 
The indices x+, X-, v, a, xt are expressed by expansions 
which are convergent for A small enough. Hence, the 
indices can be computed in principle with arbitrary pre- 
cision by an explicit computation of the first orders and 
a rigorous bound for the rest; moreover, in this way one 
can prove that the indices depend on A and on all de- 
tails of the model. On the other hand, the complexity of 
the expansions makes essentially impossible to prove the 
universal relations directly from them. 

Another important class of models whose correlations 
can be analyzed by similar methods are models of inter- 
acting fermions on a ID lattice or quantum spin chains; 



they are all described by the Hamiltonian H = 



L-l 



x=l 

^1 11 

+hY^{p^--) + X ^ v{x - y){p^ - -){py - -) 



x=l 



where are the fermion creation or annihilation op- 
erators and px = a^a^. By using the Jordan- Wigner 
transformation, the Hamiltonian of the Heisenberg quan- 
tum spin chains can be written in this way, if Ji -|- J2 = 
2, u = (Ji — 0/2)72 and J3 = —A; in particular, if 
v{x — y) = (5|a,_y|^i/2 and h = 0, we have the XYZ 
model. Let us define x = {x,xo), = e^^°Oa,e~^^° 



and, if A = Oxi...Ox„, (A) = limL^oo T be- 

ing the time order product. If u = 0, it was shown in 
[7] that, for A small enough, if T denotes the truncated 
expectation, 

^ 1 -I- Q(A) 1 -I- Q(A) 

cos[2ppx) ^^^^^^ + ivsXor]-+ ^ 2n^[x^ + {ysXa^] ^ ' 

where S^^ = p^- 1/2, pp = cos'^h + A) + 0(A) is the 
Fermi momentum (if h = 0, pp = tt/2 by symmetry ) 
and Vs = sinpp + 0(A) is the Fermi (or sound) velocity, 
which is modified by the interaction, since, contrary to 
the previous Ising case, there is no symmetry between 
space and time. Finally a;+ is a critical index, expressed 
by a convergent expansion; it depends on all details of the 
model, as it is apparent from the explicit computation of 
its first order contribution, which gives a;+ = 1 — aiA -|- 
0(A2), with ai = [v{0) - v {2p p)] /{n sin pp). When Ji ^ 
J2, that is u ^ 0, {Sx'' S^'')j, decays exponentially, with 



correlation length ^ ~ | Ji 



with z/ = 1 + aiA 



0(A ), fli being the same constant as before. If Ji 



J2, {a^apT ^ I 



-l-r, 



V 



0(A2) > 0, and the 



correlation of the Cooper pair operator p^ 
a~a~,, x' = {x + l,xo), decays as |x — y| 
l + aiX + 0(X^). 

If M = and jx = (2wF)~^[ax-i-i^x 



-2a;- 



a+a+, + 
, X- 



■ 



aZ a. 



x "'x-|-(l,0) 



], the 



following WI, for k, k-|-p close to p"^ = {uipp, 0), w = ±, 
are true: 

-iPoOp^ia^+p) + (jpvNippa^a^^p) ~ BG (7) 



k+p"k+p/ 



B 



with G = G(k, k -h p) = [( 
B = 1 + 0{X) and vj,vn ~ Vs{l + 0(A)); in particular 
vn/vj = 1 + 2aiX + O(A^), with ai the constant defined 
above, after (6). When A = 0, the continuity equations 
for px and jx imply WI similar to (7) with B = I and 
vn = vj = Vs; the interaction has the effect that the 
normalization B is not 1 ( [H, p^] = but [H, j^] ^ 0) 
and two different velocities, the charge vj and the cur- 
rent velocity vn, appear. The presence of the lattice. 
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which breaks the Lorentz symmetry of the continuum 
hmit, causes the presence of three distinct velocities, vn, 
vj. Vg. Finally, we recall that the susceptibility is de- 
fined as K = linip^o 5^(0. p), where il(po.p) is the Fourier 
transform of {S^^^S^q^} rp', Kp ^ is the compressibility if p 
is the fermionic density. Our results are contained in the 
following Theorem. 

Theorem Given the models with hamiltonian (1), (5), 
at small coupling all the indices defined above can be 
uniquely expressed in terms of one of them: 

x-=x~'^ , a = 2(l-a;+)(2-a;+)"^ (8) 
u-'^ = 2-x+ , jy-i = 2 - x^^ (9) 
2r] = x++x^^-2 (10) 

XT = {2-x+){2-x^Y' (11) 

Moreover, in the model (5) the velocities appearing in (7) 
verify VnVj = Vg and vj = sinpj?, while the susceptibility 
K verifies 

K = Xj^{'KVs)~^ ■ (12) 

(11) is a new relation for the Ashkin- Teller model, never 
proposed before; the first relation in (8) was conjectured 
in [10] and (9), (10) in [12, 17]. (12) is part of the Haldanc 
Luttinger liquid conjecture [13] for fermionic systems or 
quantum spin chains. Some of the above relations were 
checked in certain solvable case: the second of (8), which 
is equivalent, by using the first of (9), to the hyper-scaling 
relation 2v ^ 2 — a^va the case of the Eight Vertex model 
[3]; (10), (12) in the case of the Luttinger model [8]; (12) 
in the XYZ spin chain [13]. The above theorem provides 
the first proof of the validity of such relations for generic 
non solvable models. Note that, in the notation of [13], 
vn = (ttk)"^ should not be confused with vn appearing 
in the WI (7); they arc coinciding only in the special 
case of the Luttinger model. Therefore vn = vn is an 
example of relation true in the Luttinger model but not 
in the presence of a lattice. 

Outline of the proof. The technical details are long 
and will appear elsewhere [5, 7]; here we just outline 
the proof. The partition function and some of the cor- 
relations of the spin model (1) can be exactly rewritten 
as sums of Grassmann integrals describing = 1 + 1 
Dirac fermions on a lattice and with quartic non lo- 
cal (but short ranged) interactions, by using the clas- 
sical representation of the Ising model in terms of Grass- 
mann integrals [14]. The Grassmann variables are writ- 
ten as Vk — X]h=-oo ^k'^ '^^^'^ ^k*^ living at momen- 
tum scale k = 0(2'*). After the integration of the fields 
'ip^^\ . . . ,ip^^~^'^\ the partition function can be written 
[18] as 

/ Pz..,.(#(^''))e^<^'(^'^'^'") (13) 
where Pz^,ij,^{dip'^-'^^) is the Gaussian Grassmann inte- 



gration with propagator g^-'^\\s.) = ■^''^^ ^ x 

/-isinfco + sinfc -f -ph-H-+ \ ^ 
Y —ph — IJ'+- —i sin fco — sin k\-\- p J 

(14) 

where X?i(k) is a smooth compact support function non- 
vanishing only for |k| < 2^ , Zh and ph are the effective 
wave function renormalization and the effective mass, 
p±± are O(k^) and non vanishing at k = (±7r, ±7r) 
(there is no fermion doubling problem); moreover, 
'^'s V'x -i-^x -t-V'x -V'x - + Rh, where Xh is the effective 
coupling and Rh, is a sum of irrelevant terms, represented 
as space-time integrals of field monomials, multiplied by 
kernels which are analytic functions of \k, k > h. Ana- 
lyticity is a very non trivial property, obtaincxi via tree 
expansions [6] and exploiting anticommutativity prop- 
erties of Grassmann variables, via Gram inequality for 
determinants (which takes into account compensations 
among different graphs of different signs at a given or- 
der). It is important to stress that (13) is exact, in the 
sense that the irrelevant terms and the lattice are fully 
kept into account (in standard RG applications they are 
instead neglected). The effective coupling Xh converges, 
as h ^ —00, to a function X-oo (analytic function of 
A) , thanks to the asymptotic vanishing of the beta func- 
tion, which is a c;onseqiicncc of Ward Identities. A simi- 
lar analysis can be repeated in the case of the fermionic 
model (5), the main (but trivial) difference being that 
(14) is replaced by a similar expression, taking into ac- 
count that xq is a continuous variable; such asymmetry 
has the effect that, contrary to what happens in the spin 
case, the velocity is renormalized by the interaction. In 
order to exploit the asymptotic symmetries of the model, 
it is convenient to introduce the following Grassmann in- 
tegral 

I P|''(#(^^))e^'"'(^^'-"') (15) 

where, if -0 = (?/>+ , ^p- ) and = V'''7o are Euclidean 
d = 1 -|- 1 spinors, P^'^{d^^-^'>) is the fermionic gaussian 
integration with propagator ^(-^^(k) = xjv(k)(7^k^)~^, 
and T/W^V-^^^)) = /dxdyz;(x - y)j^(x).7V,(y), with 
jp(x) = tjjx'Ji^iijJx and w(x — y) a short range symmet- 
ric interaction. A multiscalc integration is now neces- 
sary also in the ultraviolet region to perform the limit 
N ^ 00, while in the integration of the infrared scales 
an expression similar to (13) is found; the effective cou- 
pling is denoted by X^- The crucial point is that it is 
possible to choose, by a fixed point argument, the values 
of Aoo (fixed c ~ sin pp in the model (1), while c = Vs in 
the model (5)) so that A_oo = A_oo- This implies that 
the critical exponents of the two models are the same, 
because the exponents are expressed by series in A_oo/c 
with universal coefficients. Of course Aqo is a convergent 
series in A depending on all details of the models (1) or 
(5). On the other hand, the contimmm Grassmann in- 
tegral (15) verifies extra Lorentz and Gauge symmetries, 
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implying exact Ward Identities when the ultraviolet cut- 
off is removed; by the transformation f/; e*"="!/'x one 
finds -ip^O'^,pV'kV'k+p)tft = 



('0k'0k)tft - ('0k+pV'l 



k+p/th 



AAr(k,p) 



(16) 



where A^r = ((5jpV'k'0k+p)t/i, with 6jp = J dk[(Xjv^(k + 

p) - l)(7pkp + 7pPm) ~ (Xjv^(k) - l)7MkM]'0k'0k+p; an 
analogous expression is obtained for the axial current 
V'7m75V'- By a multiscale analysis it can be proved that 
lim^v^oo Ajv(k,p) = 



-?:Tw(p)p^(j^,pV'k,a.V'k+p,w) 



th 



Aoo/(4^c) (17) 



A similar expression holds for the chiral WI; the fact 
that Aat (k, p) is not vanishing in the limit — s- oo is a 
manifestation of a quantum anomaly. The anomaly co- 
efficient T is linear in Aoo; this is the non-perturbative 
analogue of the anomaly non renormalization in QED in 
4D. Such crucial property depends on our assumption 
about the interaction in (15); it would not be true, for 
instance, if we replace w(x — y) with a delta function [4]. 
By combining the WI with the Schwinger-Dyson equa- 
tions, one gets some equations for the correlations, from 
which the indices can be computed as functions of r. 
One can find, for example, that x+ = (1 — t)(1 + r)^^, 
X- = (1 -f-T)(l — r)^^, so that x+X- = 1; the other rela- 
tions follow by similar arguments. Note that the_ indices 
we consider have a simple expression in terms of Aqo , but 
Aoo is of course rather complex and model dependent as 
a function of A. 

A similar RG analysis can be repeated for the 
model (5); it turns out that the vertex functions in 
the first line of (7) are asymptotically coinciding with 

Z^'HfJtJu+pJtH and iZ^'HfJtJu+^Jth^ with 



i-(3) /z(3) ^i + aiX + 0(A2); therefore, by using the WI 
for the model (15), we derive (7) with B = ^(3)^-^(1 - 
r)-i = 1, B = Z(^)Z-\1 + t)-\ vn = VsZ(^yz(^\ 
vj = WsZ^'^^/Z^'^); on the other hand, the equations of 
motion related to the lattice Hamiltonian impose the 
constraints vj = sin pp and B = 1. Finally a WI for 
the density correlation can be also derived; if -Da;(p) = 
—ipo + LVVsP and (l{p) = {ppPp) =, we get 



(Z(3))2 



47r«,Z2 1 -t2 



2- 



D_ip) D+(p) 



£>+(p) £>_(p) 



which implies (12), by using k = limp^o f^(0,p)- 

In conclusion, we have established for the first time the 
validity of a number of universal relations among critical 
exponents and other quantities in a wide class of gener- 
ally non solvable lattice models. They are true in special 
continuum solvable models and we have proven that the 
lattice symmetry breaking effects produce different veloc- 
ities in the model (5) and change the critical exponents, 
but do not destroy the validity of several universal rela- 
tions (on the other hand, not all the relations valid in 
the solvable models are generically true, like the relation 
Vn = Vn ot the relations for the dynamic exponents [15]). 
Some of the universal relations are used for the analysis 
of experiments in carbon nanotubes or spin chains, but 
we believe that their interest goes much beyond this, as 
they provide one of the very cases in which the universal- 
ity principle, a general belief in Statistical Physics and 
beyond, can be rigorously verified. Extensions of our 
methods will hopefully allow to prove universal relations 
in an even wider class of models, as well as other relations 
among spin or dynamic exponents. 
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